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A further comparison of graphene and thin metal 
layers for plasmonics 


Xiaoyong He,** Pingqi Gao” and Wangzhou Shi® 


Which one is much more suitable for plasmonic materials, graphene or metal? To address this problem 
well, the plasmonic properties of thin metal sheets at different thicknesses have been investigated and 
compared with a graphene layer. As demonstration examples, the propagation properties of insulator— 
metal—insulator and metamaterials (MMs) structures are also shown. The results manifest that the plasmo- 
nic properties of the graphene layer are comparable to that of thin metal sheets with the thickness of tens 
of nanometers. For the graphene MMs structure, by using the periodic stack structure in the active region, 
the resonant transmission strength significantly improves. At the optimum period number, 3-5 periods of 
graphene/SiOz, the graphene MMs structure manifests good frequency and amplitude tunable properties 
simultaneously, and the resonant strength is also strong with large values of the Q-factor. Therefore, 
graphene is a good tunable plasmonic material. The results are very helpful to develop novel graphene 


plasmonic devices, such as modulators, antenna and filters. 


1. Introduction 


Surface plasmons (SPs) are two-dimensional (2D) surface 
waves propagating along the boundary between a metal and 
dielectric medium, offering the merits of controlling the light 
at sub-wavelength scales.” They exploit the mass inertia of 
electrons to create propagating charge density waves at the 
metal surface and exhibit sub-wavelength confinement of the 
electric field in the direction normal to a conductor-dielectric 
interface.** As a novel branch of electromagnetic structures, 
plasmonics has emerged as an important scientific field, 
promising bright prospects in the aspects of waveguide 
devices, solar cells, and spectroscopy.” ** Considerable efforts 
have been devoted to develop various plasmonic systems, i.e. 
the metal-dielectrics-metal structure, the groove structure," 
hybrid plasmonic structures,*° and novel metamaterials (MMs) 
devices. To find suitable material systems is of vital impor- 
tance for the further development of plasmonic systems. As 
the fundamental plasmonic materials, noble metals (Au and 
Ag) exhibit the disadvantages of short wavelength, not very 
good confinement and inconvenient tunability. Thus, other 
kinds of optional materials have been proposed, such as 
graphene, semiconductors,” superconductors’ etc. 
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Graphene, a flat monolayer of carbon atoms packed into a 
honeycomb crystal lattice,” is a fascinating material and 
shows the advantages of high carrier mobility, flexible tunabil- 
ity without varying physical geometrical parameters, and 
strong mode confinement.” Compared with its metal 
cousins, graphene is a good complementary material for the 
exploration of SPs and MMs devices.” *? Much experimental 
and theoretical research work has been carried out in this 
aspect. For instance, Xia’s group experimentally investigated 
the damping mechanisms of GSPs in the mid-infrared regime, 
indicating that the surface polar phonons in the SiO, substrate 
resulted in strong plasmonic dispersion, very useful in design- 
ing novel graphene plasmonic devices.*” By depositing metal 
sub-wavelength ring apertures on the graphene-SiO,-Si layers, 
Gao et al. suggested a high-speed modulator, indicating that 
the amplitude modulation of THz waves can reach about 
50%.°* With the help of pumped green light, a terahertz active 
diode has been experimentally demonstrated by using a hybrid 
graphene-Si structure,** which permits THz waves to pass 
through at a positive biased voltage and blocks the waves 
under a small negative voltage. The corresponding trans- 
mission modulation depth reaches more than 80%, indicating 
good applications in modulators and switches. 

Probably, one of the most important challenges to the 
applicability of plasmonic devices is to exhibit good confine- 
ment and small dissipative losses (long propagation length) 
simultaneously. To find a suitable material is an important 
step to solve this problem. As the strong candidate, graphene 
or metal film, which one is more suitable for plasmonic 
materials? Some researchers claim that metal is superior to 
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graphene membrane for plasmonic devices.** But if thin metal 
sheets (tens of nm) are adopted as the composed materials, 
large errors between the experimental and simulation results 
occur on the basis of bulk metal parameters.** Metallic nano- 
scale sheets have also become a good platform to investigate 
plasmonic devices.** Because of the surface scattering and 
backscattering, the plasmonic properties of a thin metal film 
closely depend on structural dimensions if its thickness is 
comparable to the electron mean free path (MFP).°°** To a 
large extent, it is unfair to compare the plasmonic properties 
of a graphene layer with those of the bulk metal database. 
Therefore, the plasmonic properties of a graphene layer and 
thin metal sheets have been theoretically investigated and 
compared. Furthermore, as demonstration examples, the 
propagation properties of insulator-metal-insulator (IMI) and 
sub-wavelength MMs structures based on graphene and met- 
allic MMs patterns have also been shown. The results manifest 
that the plasmonic properties of the graphene layer are com- 
parable to those of the thin metal membrane if the thickness 
of the metal film is in the range of tens of nanometers. 


2. Research methods 


Fig. 1(a) and (b) show the side view of mono and multi-periods 
of graphene (metallic) pattern MMs structures. The substrate 
is made from polyimide with a thickness of 2 um. Fig. 1(c) and 
(d) display the top views of electric split-ring resonators and 
circular ring shaped unit cell structures, respectively. The 
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incident waves are normally transmitted through the multi- 
layer structure along the z direction. 

For the bulk model, the frequency-dependent permittivity 
of the metal in the THz regime is as follows:*° 


Wp” 
olw + iwt) 


elw) = £% (1) 
For thin metal sheets, owing to the surface scattering, the 
conductivity of a thin metal film strongly depends on the 
thickness. According to the Fuchs-Sondheimer (FS) theory and 
simplifying the Boltzmann equation in a scalar relaxation time 
approximation, the ratio between the resistivity p and the 
intrinsic film resistivity p; (Łe., the resistivity the film would 
possess if it is infinitely thick) is given by:*° ** 
Pi i 
p 


= 1.0 


(3/8x) 4 (1.0/t?-1.0/t)e“dt (2) 


(3 /2) | 
1 
where x (=d/l;) is the ratio between the film thickness d over 
the intrinsic electron mean free path li. 
In addition, according to the Drude-Smith (DS) model, 
taking into account the back-scattering, the conductivity of 
thin metal films is expressed as:*”*° 


E00@p?T C; 
1.0 S k 2 (1-— = (3) 


Besides the common Drude term, the second part in eqn (3) 
gives the persistence of electron’s velocity after a number of 
scattering events, the parameter j denotes the scattering 
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(a), (b) The side view of the MMs structure. (a) The graphene unit cell structure is deposited on SiO2/Si layers, (b) the multilayer metal—dielec- 


trics structure, the thickness of the dielectric layer is 30 nm, and the doped Si layer is used to apply the gate voltage. (c), (d) The top views of geome- 
try for the MMs unit cell structures. (c) The electric split ring resonator (eSRR) structure, lx = ly = 48 um; (d) the circular ring structure, Dı = 40 um, D2 


= 36 um. The periodic length along the x and y directions are both 60 um. 
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events, c; describes the persistence of velocity, i.e. the fraction 
of electron’s initial velocity that is retained after the j® 
collision. 

As graphene is a typical 2D material, its surface conductivity 


can be described by o,:** 


Og(@, Me, T, T) = Ointer + Ointra = Gel x 
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in which fa(e) is the Fermi-Dirac distribution, ¢ is the energy of 
the incident wave, kg is Boltzmann’s constant. The first and 
second parts denote intra-band and inter-band 
contributions. 

As a typical example, the sub-wavelength insulator-metal- 
insulator (IMI) structure has been adopted to compare the 
graphene layer and thin metal layer. For thin metal layers, 


the 


the dispersion relationship for the IMI structure can be 


~ y (24) raf. (5) 


in which £m and eg are the dielectric constants of the thin 
metal layer and surrounding dielectric constants, A is the 
thickness of the thin metal layer. For the graphene layer, the 


written as:*” 


dispersion relationship can be expressed as:**”° 
& & 6 
1474 5-8 = 9, (6) 
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where œ<; is the dielectric constant of the surrounding 
materials, # is the propagation constant, y;” = J” — ko"e;, and og 
is the graphene conductivity. 

Furthermore, for the multilayer structure, a numerical 
transfer matrix method (TMM) is applied to study the 
propagation properties. The relationship between the magni- 
tude of the magnetic field of incident light at the first layer 
and that of the reflected waves at the last layer are shown as 
follows:*° 
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Mi; is the characteristic matrices in the i® layer, d; is the thick- 
ness of the i” layer. 

The surface plasmon wavelength 4,,, propagation length Lsp 
and decay length Dsp have been defined as:"” 


Lspp = 1/|Im(A)| (9) 


Dspp = 1/|Re( VA =k?) 


3. Results and discussion 


Fig. 2 shows the permittivity of the monolayer graphene mem- 
brane and metal sheets with different thicknesses. The permit- 
tivity of the thin metal layer was obtained from the Drude, 
Fuchs-Sondheimer and Drude-Smith models. The solid black 
and orange color lines are for the permittivity of the bulk 
metal layer (Au) and graphene; the other solid color lines are 
the results obtained from the FS theory; the dashed lines show 
the results from the DS model. Au nano-sheets are used 
because of their good conductivity and chemical inertness. 

As shown in Fig. 2, for the real part of permittivity of 
metallic thin films, the influence of thickness on plasmonic 
properties is significant. If the thickness is very small, e.g. 
5 nm and 10 nm, their metal permittivity is much smaller 
than that of bulk metal films. As the thickness increases, the 
value of metal permittivity increases drastically, which means 
that the metal layer denotes much better conductive pro- 
perties. If the metal sheet thickness is larger than 100 nm, the 
distinctions between different models become smaller. If the 
thickness reaches more than 200 nm, the metal layers can be 
regarded as bulk materials. For instance, at 1 THz frequency, if 
the metal layer thicknesses are 10 nm, 20 nm, and 50 nm, the 
ratios of their permittivity to bulk value are 0.36, 0.52, and 
0.73, respectively. For the results of the graphene layer, the 
value of Re(e,) is larger than that of metal nano-sheets if the 
operation frequency is low, which means that the graphene 
layer is very a good plasmonic material in the THz regime. 
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Fig. 2 Comparison of the dielectric constant of graphene layer (Eş = 0.3 
eV) with the metal films at different thicknesses. The solid and dashed 
lines are for the FS and DS models, respectively. The thicknesses of the 
thin metal layers are 5 nm, 10 nm, and 20 nm, respectively. 
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In addition, the imaginary parts of the graphene layer can be 
found in Fig. 2(b), which manifest that the values of Im(e,) 
(orange line) are roughly equal to those of thick metal sheets 
(black solid lines) if the operating frequency is larger than 
1 THz. Thus, overall we can infer that the dielectric properties 
of the graphene layer are comparable to those of the thin 
metal layer in the THz and mid-infrared spectral regime. 

To give a good understanding and comparison of the plas- 
monic properties of thin metal layer and graphene, as an 
example, the dispersion relationship based on the IMI struc- 
ture has been shown in Fig. 3. The Fermi level of the graphene 
layer is 0.3 eV. The thicknesses of the graphene layer and thin 
metal layer are both 0.34 nm. The dielectric constants of thin 
metal layers are obtained from the FS model. The dispersion 
properties for graphene and thin metal layers have been deter- 
mined by using eqn (5) and (6), respectively. As shown in 
Fig. 3(a), the real part of the effective index of the thin metal 
layer, Nefm, is larger than that of the graphene layer, Nerf, 
which indicates that the thin metal layer shows better confine- 
ment for the IMI structure. Nevertheless, the losses of the thin 
metal layer are also much higher, the values of Im(“efrm) are 
at least two orders larger than those of graphene IM(neffg). 
This phenomenon is closely associated with the dielectric pro- 
perties of thin metal layers and graphene. For instance, if the 
operation frequencies are 1 THz, 10 THz, and 30 THz, the 
dielectric constants of the graphene layer are —1.10 x 10° + 
6.59 x 10°i, —1.46 x 10* + 9.18 x 1077, and —1.56 x 10° + 3.54 x 
10*i, respectively. The values for the thin metal layer are —1.90 
x 107 + 4.18 x 1077, —1.05 x 10° + 2.32 x 1077, and —1.12 x 107 + 
8.95i, respectively. Furthermore, the ratio of Re(€,)/Im(em) is 
also a common tool to measure the “metallic” properties of 


ae Re(k === Im(k,) 
=—Re(k = Im(k 
elk == IME am 


— =e ~ 


p 40 60 


— 
2. 


= Graphene 
—— Metal 


Figure of merits 
2 


— = 
2 Q 
D p 


20 
Frequency (THz) 


40 60 


Fig. 3. (a) The effective indices and (b) the figures of merit of graphene 
and thin metal film based on the IMI structures. The Fermi level of the 
graphene layer is 0.3 eV. The thicknesses of thin metal layer and gra- 
phene are both 0.34 nm. 


plasmonic materials. At the frequency of 30 THz, the values of 
Re(ém)/Im(ém) for graphene and thin metal layer are 44.2 and 
12.5, respectively. Obviously, the values of dielectric constant 
and Re(é)/Im(em) of the graphene layer are much larger than 
those of the thin metal layer, which implies that the graphene 
layer is a better “plasmonic” material. In addition, the figure- 
of-merit (FOM) for confinement-to-propagation can also be 
found in Fig. 3(b), which has been defined as Re(/)/Im(f), 2 is 
the propagation constant. For the thin metal layer, as the fre- 
quency increases, the real part of ne increases, resulting from 
the smaller values of the dielectric constant at larger fre- 
quency, while the imaginary part of ne changes little. There- 
fore, the FOM of the thin metal layer increases with frequency 
gradually. But if the operation frequency is high, larger than 
80 THz, the values of FOM decrease due to the significant 
increase in the loss of metal, as shown in Fig. 3(a). While for 
the case of the graphene layer, as the operation frequency 
increases, the real and imaginary parts of me both increase, 
and the values of IM(neffg) increase much faster at low fre- 
quency, leading to a dip of FOM at low frequency. Compared 
to the values of the thin metal layer, the FOM of the graphene 
layer is much larger because of the smaller losses of the 
graphene layer. For example, at the frequency of 1 THz, 
10 THz, and 30 THz, the FOM values of the graphene layer are 
49.2, 20.2, and 45.2, and the FOM values of the thin metal 
layer are 0.469, 4.52, and 12.5, respectively. Therefore, from an 
overall viewpoint, a graphene layer is a much better plasmonic 
material compared to the thin metal layer. 

The plasmonic properties of multilayer structures for 
graphene and metal layers with different thicknesses are 
shown in Fig. 4, which have been obtained by using the trans- 
fer matrix method. The inset in Fig. 4(a) shows the multilayer 
waveguide structure. The refractive indices of SiO2, Al,O3 and 
Si are 1.98, 2.7, and 3.42, respectively. The permittivity of the 
metal films was obtained from the FS and DS models, respec- 
tively. As the thickness of dielectric filling materials increases, 
the propagation length increases, resulting from the decrease 
of losses for the waveguide structure. In addition, as the metal 
film thickness increases, the propagation length increases. 
This is explained by the fact that with the decrease of 
thickness, the metal permittivity reduces, displaying poor con- 
ductive properties, leading to larger losses and smaller propa- 
gation length. For example, at the frequency of 3 THz, the 
thickness of the dielectric stripe is 30 pm, while the metal 
layer thicknesses are 10 nm, 20 nm, and 50 nm, the permittiv- 
ity values of the metal film are —8.48 x 10* + 6.20 x 10%, 
—1.21 x 10° + 8.82 x 107, and —1.69 x 10° + 1.24 x 10°7. The 
corresponding propagation lengths are 2.52 x 10° um, 4.83 x 
10° um, and 1.08 x 10* um, respectively. If the metal film thick- 
ness is 50 nm, the results of metal films approximate to those 
of the bulk metal, which is also in accordance with the permit- 
tivity of thin metal layers, as shown in Fig. 2. 

For the case of a graphene layer, the propagation length is 
in the scope of thousands of micrometers, slightly smaller 
than that of a metal film. As the Fermi level increases, the 
carrier concentration and permittivity of the graphene layer 
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Fig. 4 Comparison of (a) the propagation length (b) decay length plas- 
monic properties of thin metals with different thicknesses and graphene. 
The Fermi levels of graphene layer are 0.5 eV and 1.0 eV. The thicknesses 
of the Au layer are 10 nm, 20 nm, 50 and 100 nm, respectively. The 
operating frequency is 3.0 THz. The inset in (a) shows the multilayer 
waveguide structure. 


increases, resulting in the enhancement of losses. The lateral 
decay length is used to measure the confinement of propa- 
gation modes, as shown in Fig. 4(b). With the increase of the 
dielectric stripe thickness, the decay length increases, indicat- 
ing poor mode confinement. The decay length decreases with 
the thickness of the thin metal layer. For the graphene layer, 
the decay length is much smaller than that of the metal film, 
which means that the graphene layer shows much better con- 
finement than those of the metal film. 

Additionally, for the practical application of plasmonic 
devices, good confinement usually results in large losses and 
small propagation length. To give an overall benchmark of the 
plasmonic properties of the above-mentioned materials, the 
values of figure-of-merit versus dielectric filling material thick- 
ness are shown in Fig. 5, which can be defined as Re(/)/Im(f). 
The structural parameters are the same as those in Fig. 4. If 
the metal film thickness is more than 50 nm, the FOM is in 
the range of several hundreds, much larger than that of the 
graphene layer. With the decrease of the metal layer thickness, 
the value of FOM decreases. If the metal sheet thickness is 
20 nm, the ratio values of the metal layer are roughly equal to 
those of the graphene layer. If the metal layer is thinner, 
<20 nm, the values of FOM for the thin metal layer are smaller 
than those of the graphene layer, indicating that the plasmo- 
nic properties of the graphene layer are better than those of 
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Fig. 5 Comparison of the figure-of-merit for metals thin layers and 
graphene versus the thickness of the dielectric filling materials. The 
Fermi levels of graphene layer are 0.5 eV and 1.0 eV. The thicknesses of 
the Au layer are 10, 20, 50, and 100 nm, respectively. 


the metal layer. Thus, we can draw a conclusion that the plas- 
monic properties of the graphene layer are comparable to 
those of thin metal layers with the thickness in the range of 
tens of nanometers. This is different from the conclusion in 
many published articles, where the authors claim that the 
metal system is much better than those of graphene based 
devices. The metal layer thickness is very small in those 
articles, in the range of tens of nanometers, but they still use 
the bulk metal database, not including the influences of metal 
layer thickness. Actually, the permittivity of the thin metal 
layer tightly depends on thickness, especially when the value 
of thickness is comparable to mean free path. When we 
discuss and compare the plasmonic properties of the metallic 
plasmonic systems, the influence of metal sheet thickness 
should be taken into consideration. 

With the designed electromagnetic properties, sub-wave- 
length MMs structures are a good application playground for 
plasmonic devices. Recently, the graphene layer has been inte- 
grated with the MMs structure to develop tunable devices. To 
compare the plasmonic properties of graphene and thin metal 
sheet further, Fig. 6 shows the transmission curves of the MMs 
structures. The simulation results were acquired by using the 
CST Microwave Studio. The unit-cell boundary conditions in 
the x-y plane and Floquet ports in the z direction have been 
adopted. By using the obtained S-parameters from the simu- 
lation, the transmission (T(q)) curves can be achieved by the 
formula, T(@) = |So;|’. A typical MMs element, electric split 
ring resonator structure is adopted as the unit cell, as shown 
in Fig. 1(c). The polarization of incident light is along the 
y direction. The solid black line is for the results obtained 
from the bulk Drude model, the metallic eSRR thickness is 
0.2 um. The transmission curves of the bulk metal shows two 
dips, i.e. the dips at the frequencies of 0.95 THz and 3.02 THz, 
which correspond with the LC and dipolar resonances, 
respectively. 
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Fig. 6 (a) The transmission of MMs structures based on the eSRR unit 
cell structure for different kinds of materials. The polarization of the 
incident wave is along the y direction. (b) The polarization is along the x 
direction. The thicknesses of the metal layer are 10 nm, 20 nm, and 
50 nm, respectively. The Fermi levels of the graphene layer are 0.5 and 
1.0 eV, respectively. 


As the sheet thickness decreases, the permittivity of the 
metal reduces, the resonant strength of transmission curve 
becomes weak, and the resonant dip moves to low frequency 
as well. For instance, under the FS model, if the metal film 
thicknesses are 10 nm, 20 nm, and 50 nm, the values of trans- 
mission dips (resonant frequencies) for the LC resonance are 
0.36 (0.82 THz), 0.17 (0.90 THz) and 0.062 (0.93 THz). For high 
frequency dipolar resonances, the values of transmission dips 
are 0.052 (2.17 THz), 0.008 (2.94 THz) and 5.03 x 10~* (2.97 
THz). The dashed lines show the numerical results from the 
DS model. Next, we discuss the impact of metal sheet thick- 
ness on the transmission curves. It can be found from Fig. 6(a) 
that if the metal film thickness is 10 nm, there is no trans- 
mission dip because the metal sheet shows very poor conduc- 
tive properties. As the metal film thickness increases, reaches 
the 20 nm, the conductive properties of the metal film 
increase, and the obvious resonant dip appears at the low and 
high frequencies. If its thickness increases further, reaches 
more than 50 nm with the metal layer showing good conduc- 
tive properties, the corresponding permittivity can reach about 
70% value of the bulk metal. The metallic MMs structure at 
this thickness manifests good resonance properties. In a word, 
if the metal film thickness is lower than 10 nm, it cannot show 
an obvious transmission dip; when its thickness reaches more 
than 50 nm, the MMs structure shows strong transmission 
dips, approximating the properties of the bulk metal layer. In 
addition, for the graphene pattern structures, the influences of 
Fermi levels on the transmission curves have also been shown. 


If the Fermi level of the graphene layer is high, ie. 1.0 eV, the 
sub-wavelength MMs structure shows an obvious dip at low fre- 
quency (LC resonance, 0.68 THz) and high frequency (dipolar 
resonance, 1.78 THz). With the decrease of the Fermi level, the 
transmission dip at lower frequency does not seem to be 
evident. 

Fig. 6(b) shows the resonant transmission spectral curves if 
the polarization is along the x direction. There is no low fre- 
quency transmission dip, different from the result that the 
polarization direction is along the y direction, but the broad 
transmission dip at high frequency still exists. Even if the 
metal film thickness is 10 nm, the transmission dip is 
obvious, the value of the transmission dip is 0.027 (1.38 THz). 
If the metal film thickness is larger than 20 nm, the trans- 
mission dips locate near 1.5 THz. This is due to the fact that 
dipolar resonance mainly depends on the structural para- 
meters, unlike the LC resonance, very sensitive to the conduc- 
tive film thickness. For the graphene pattern MMs structure, 
the transmission resonance becomes weaker, and the resonant 
dip shifts to low frequency with the decrease of the Fermi 
level. 

The 2D field plot is a good way to understand the resonant 
mechanisms of MMs structures. The magnetic component 
along the z direction and surface current density can be seen 
in Fig. 7. The polarization direction is along the y direction. 
The materials are the bulk metal, metallic thin sheet with the 
thickness of 20 nm, and the graphene layer with the Fermi 
level of 1.0 eV. Fig. 7(a)-(d), and (e), (f) are the results for the 
bulk metal, metal layer with the thickness of 20 nm, and the 
graphene membrane, respectively. The resonant frequencies 
are 0.95 THz, 0.90 THz and 0.68 THz, respectively, which is 
accordance with the LC resonant frequency, as shown in 
Fig. 6(a). In Fig. 7, red indicates high intensity, and blue is for 
low intensity. The LC resonance results from the fact that the 
incident THz waves interacts with the magnetic modes of the 
MMs structure. The results are similar to Fig. 2(e) and (g) in 
ref. 33. It can be found from Fig. 7(a), (c) and (e) that with the 
decrease of the membrane thickness, the value of the mag- 
netic field decreases, i.e. the maximum values of H, are 1.95 x 
107, 1.03 x 10* and 5.59 x 10° Vm’, respectively. This denotes 
that the transmission resonant strength decreases with the 
sheet thickness. Fig. 7(b), (d) and 7f) show the corresponding 
results of the surface current density, which manifests that the 
incident THz fields drive the circulating surface currents in 
the inductive loops, inducing the charge accumulation at the 
capacitive split gaps. The surface current figures are very 
similar to the results shown in Fig. 1(c) in ref. 33, which con- 
firms that the resonance at the low frequency (0.76 THz) is LC 
resonance. As the thickness of the metal film decreases, 
the relatively absolute values of the surface current density 
decreases, and the values are 218.1, 208, and 191.3 A m°, 
respectively. The relatively low values of the surface current 
density mean small resonant strength, which corresponds to 
the results shown in Fig. 6(a). 

In accordance with the high frequency of the dipolar reso- 
nance shown in Fig. 6(a), Fig. 8 shows the magnetic com- 
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Fig. 7 Shows the H; components of and the surface current density of the eSRR-shaped MMs structure. The polarization direction of the incident 
light is the along y direction. (a), (b) The bulk metal, (c), (d) are the result for the metal film with the thickness of 20 nm, and (e), (f) and are the gra- 
phene eSRRs with the Fermi level 1.0 eV, respectively. The resonant frequencies are 0.95 THz, 0.90 THz, 0.68 THz, respectively. 


(a) (c) 


Fig. 8 


(e) 


(a)-(f) Show the H, components of and surface current density of the eSRR-shaped MMs structure. The polarization direction of the incident 


light is along the y direction. (a), (b) The bulk metal, (c), (d) are the results for the metal film with the thickness of 20 nm, and (e), (f) are the graphene 
eSRRs, respectively. The resonant frequencies are 3.03 THz, 2.94 THz and 1.78 THz, respectively. 


ponents along the z direction and the surface current density 
for different kinds of materials if the polarization direction is 
along the y direction. The materials are the bulk metal layer, 
metallic film with a thickness of 20 nm, and graphene pat- 
terns. The permittivity of the metal layer was obtained from 
the FS model. Fig. 8(a)-(d), and (e), (f) show the results of the 
bulk metal, metal layer with the thickness of 20 nm. Fig. 8(a), 
(c) and (e) display that with the decrease of metal film thick- 


ness, the value of magnetic fields decreases. For instance, the 
maximum values of H, are 9.28 x 10°, 8.1 x 10° and 6.46 x 10° 
V m™t, respectively. Fig. 8(b), (d) and (f) display the surface 
current density, which indicates that the incident THz field 
drives the circulating surface currents moving along the 
surface. The above results confirm that the resonances at high 
frequency are dipolar resonance. Furthermore, with the 


decrease of metal film thickness, the value of surface current 
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density decreases. The relatively low value of surface current 
density means small transmission resonant strength, which 
agrees with the results shown in Fig. 6(a). 

Fig. 9(a) shows the transmission curves of sub-wavelength 
MMs structures based on the metal circular ring unit cell. The 
outer and inner diameters of the circular ring structures are 
40 pm and 36 um, respectively. The solid black line is for the 
bulk metal layer. The color lines are for the metallic sheet 
structure at different period numbers, and the thickness of 
each metallic sheet is 20 nm. As shown in Fig. 9, even if only 
one circular ring layer is used, the metallic MMs structure can 
manifest good resonance properties. Its value of transmission 
curve dip is 0.08 close to the frequency of 1.81 THz. This can 
be explained by the relatively large permittivity of the metal 
membrane at the thickness of 20 nm. For example, near the 
resonance frequency, the permittivities of the bulk metal and 
thin metal sheet (20 nm) are —3.82 x 10° + 4.46 x 10°, —1.99 x 
10° + 2.32 x 10°, the ratio of the permittivity of the thin metal 
sheet layer (20 nm) to that of the bulk metal is 0.52. As the 
period number increases, the resonant transmission becomes 
strong, and the resonant dip shifts to higher frequency. If 
for five periods Au-SiO, thin sheets are used, the metallic 
MMs structure shows similar properties to the bulk metal 
structures. 

Because of the small thickness, the conductive properties 
of the graphene layer are poor, limiting its improvement of 
modulation depth. This problem may be solved by using the 
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Fig. 9 The transmission of MMs structures based on the circular ring 
shaped unit cell structure for different layers. (a) The metallic MMs struc- 
ture with the layer number of 1, 2, 3, and 5. The metal layer thickness is 
20 nm. (b) The graphene MMs structure at the Fermi levels of 0.1 eV and 
1.0 eV. The layer number of MMs structure is 1, 2, 4, 8, 16 and 20 layers. 
The inset in (a) shows the circular ring unit cell structure. The influence 
of the period number of the graphene patterns on modulation depth 
can be found in the inset of (b). 


stack periodic structures, composed of isolated monolayer gra- 
phene separated with a thin dielectric film (e.g. SiO), as 
shown in Fig. 1(b). The transmission curves for the graphene 
multilayer MMs structure are shown in Fig. 9(b). The Fermi 
levels of the graphene layer are 0.1 eV and 1.0 eV, respectively. 
The period numbers of graphene-SiO, multilayer structure are 
1, 2, 4, 8, 16, and 20, respectively. The transmission resonant 
curve shows an obvious blue shift with the increase of the 
Fermi level. For the monolayer graphene, if the Fermi level 
changes in the range of 0.1-1.0 eV, the values of transmission 
resonance dip change in the scope of 0.28-0.82, and the reso- 
nant frequency modulates in the range of 0.40-1.08 THz. The 
modulation depths of transmission (Toa) and frequency 
(fmoa) are about 65.8% and 62.9%, which have been defined 
as Tmoa = AT/Tmax and fmoa = Aflfmax, respectively. In the THz 
regime, the intra-band conduction dominates, and the carrier 
concentration of the graphene layer increases significantly 
with the increase of the Fermi level, so the value of kinetic 
inductance of the graphene unit cell ZL, and total inductance 
Lsum decreases significantly. Consequently, the resonant dip 
shifts to high frequency and can be tunable in a wide range by 
using the gate voltage. Additionally, the effects of the layer 
number of the graphene patterns in the active region are 
shown in Fig. 9(b). As the period number of graphene patterns 
layers increases, the active area of MMs structure increases, 
and the values of kinetic inductance decreases, leading to the 
resonant frequency of transmission curve moving to high fre- 
quency. For instance, if the period number of graphene 
pattern layers is one and ten, the tunable frequency ranges are 
0.40-1.08 THz and 0.92-1.63 THz, respectively. The value of 
transmission decreases with the increase of period number. If 
the period number of graphene patterns reaches 20, the net 
transmission achieved is close to the frequency of 1.81 THz. 
The inset in Fig. 9(b) shows the frequency modulation fmoa 
and amplitude modulation Tmoa versus the period number of 
graphene pattern layers in the active region. As the period 
number increases, the fmoa decreases, while the modulation 
depth of transmission increases and reaches a relatively satu- 
rated value of about 95%. From the viewpoint of practical 
application, it is expected that the tunable graphene based 
devices manifest good frequency and amplitude modulation 
simultaneously. We define the parameter I as: I = fmoa X Tmoa- 
The values of I are 0.409, 0.452, 0.436, 0.367, and 0.258, 
respectively, when the period numbers are 1, 2, 4, 8, and 20. 
This means that as the period number of graphene layers 
increases, the stack MMs structure shows much better tunable 
properties, and the values of frequency and amplitude modu- 
lation depth are both large. However, if the period number of 
graphene pattern layers increases further, larger than 4, 
though the values of Tmoa increases, finoa and the tunable 
range of resonant frequency and amplitude decreases. 

The Q-factor is an important means to measure the reso- 
nant quality of the spectral curves, which can be defined as: 
Q = W&o/A@gwum, Where wo is the resonant frequency of the 
transmission curve dip, FWHM means the full width half 
maximum bandwidth. The values of the Q-factor for the 
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metallic thin sheets and graphene patterns at different period 
numbers are shown in Tables 1 and 2, respectively. For the Au 
circular ring structure, as the period number increases, the 
resonant frequency of the transmission dip shifts to high fre- 
quency, the spectral curve FWHM decreases, leading to a 
slight increase in the values of the Q-factor. If the period 
number is five, the resonant frequency and Q-factor are 1.85 
THz and 1.75, approximating the results of the bulk Drude 
metal layer with the values of 1.87 THz and 1.90, respectively. 
As a whole, the metallic thin sheet structures demonstrate 
good resonant properties, even if one-period is used. 

The results of graphene patterns composed of a stack struc- 
ture are shown in Table 2. If the Fermi level of the graphene 
layer is small, e.g. 0.1 eV, the circular ring MMs structure 
does not display very good resonant properties because of the 
poor conductivity of the graphene layer. The corresponding 
Q-factor of the spectral curve is merely about 0.765, if one 
period of graphene patterns is adopted, indicating the weak 
resonance strength, represented as the black dashed lines 
shown in Fig. 9(b). As the period number of graphene patterns 
increases, the resonant frequency and FWHM increases, the 
resonant strength of the spectral curves also become strong, 
leading to an increase in the Q-factor. While for the results of 
the large Fermi level of the graphene layer, e.g. 1.0 eV, the cir- 
cular structure indicates good resonance properties even if 
only one period of graphene patterns is employed, and the 
Q-factor of the transmission curve can reach more than 1.45. 
As the period number increases, the resonant frequency, and 
FWHM, and the Q-factor of the spectral curve increases. On 
the contrary, while the period number increases further, larger 
than 10, the Q-factor decreases. The reasons may be due to the 
fact that if the period number of the graphene layer is larger 
enough, the losses also improved drastically, resulting in a 
decrease in the values of Q-factors and FWHM, as shown in 
Table 2. The tunable properties of graphene supported MMs 


Table 1 The comparison of Q-factor of the different period layers of 
metallic thin layers 


Layer number Resonance frequency FWHM Q-factor 
Drude 1.87 0.99 1.90 
1 1.82 1.22 1.48 
2 1.83 1.13 1.62 
3 1.84 1.10 1.67 
5 1.85 1.06 1:75 


structure improves significantly by using the stack multilayer 
structure. If the period number of graphene patterns in the 
active region is too large, the fabrication process becomes 
difficult. Furthermore, the dissipative effects of the graphene 
patterns also become drastic, leading to the modulation pro- 
perties (fnoa and Tmoa), and the Q-factor of the spectral curves 
decreases as well. Considering the results shown in Fig. 9, the 
optimum period number of graphene pattern layers is about 
3-5, in which both frequency and amplitude modulation are 
achieved, and the Q-factor of the resonant spectral curve is 
also relatively large. 


4. Conclusion 


In summary, to find suitable materials for plasmonic devices, 
by using several kinds of methods, i.e. Drude, DS, and FS 
models, the plasmonic properties of thin metal sheets at 
different thicknesses have been investigated and compared 
with the graphene layer. As demonstration examples, the 
propagation properties of IMI and sub-wavelength MMs struc- 
tures made of graphene and thin metal layers have also been 
given and compared. The results demonstrate that the plasmo- 
nic properties of the graphene layer are comparable to those of 
thin metal sheets with the thickness of tens of nanometers. 
For the example of sub-wavelength MMs resonance structure, 
if the metal layer thickness is small (<10 nm), it only manifests 
weak metallic properties and cannot excite an obvious trans- 
mission resonance. If the metal sheet thickness is about 
20 nm, it shows relatively good metallic properties (its permit- 
tivity can reach more than half the value of the bulk metal) 
and an obvious resonant dip is observed. On condition that 
the metal layer thickness is thick enough (>50 nm), it mani- 
fests relatively much better metallic properties than the gra- 
phene layer. For the case of the graphene layer, though its 
thickness is small, the influence on the propagation properties 
is well evident. Via varying the Fermi level of the graphene 
layer, the resonant transmission curves can be tunable in a 
broad range, e.g. the modulation depth of transmission curves 
can achieve more than 60% if the Fermi level changes in the 
scope of 0.1-1.0 eV. Furthermore, by using multiple periods of 
the graphene pattern structure, the transmission resonance 
and tunable properties of the MMs structures improve signifi- 
cantly, and the net transmission is achieved if the period 
number of graphene patterns is about 20. The optimum 


Table 2 The comparison of Q-factor of the different period layers of graphene patterns in the active region 


Layer Resonance frequency FWHM Q-Factor Resonance frequency FWHM Q-Factor 
number (0.1 eV) (0.1 eV) (0.1 eV) (1.0 eV) (1.0 eV) (1.0 eV) 
1 0.41 0.528 0.765 1.08 0.746 1.45 

2 0.58 0.595 0.968 1.32 0.788 1.67 

4 0.77 0.658 1.173 1.51 0.840 1.79 

8 0.99 0.703 1.406 1.63 0.886 1.84 

16 1.25 0.798 1.560 1.78 0.969 1.84 

20 1.32 0.819 1.614 1.81 0.991 1.83 
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period number of graphene pattern layers lies between 3-5, 
where the MMs structure manifests good frequency and ampli- 
tude modulation properties simultaneously, and the resonance 
strength is also strong with larger values of the Q-factor. There- 
fore, graphene is a good tunable alternative plasmonic 
material. The results are very helpful to develop novel gra- 
phene plasmonic devices (e.g. modulators, antenna and filters) 
and find potential applications in the fields of biomedical 
sensing and wireless communications. 
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